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Abstract 

Four series of tensile loading-unloading tests are performed on isotactic polypro- 
pylene in the sub-yield domain of deformations at room temperature. In the first 
series, injection- molded specimens are used as produced, whereas in the other series 
the samples are annealed for 24 h at 120, 140 and 160 °C, which covers the low- 
temperature region and an initial part of the high-temperature region of annealing 
temperatures. A constitutive model is developed for the elastoplastic behavior of 
a semicrystalline polymer. The stress-strain relations are determined by five ad- 
justable parameters that are found by fitting the experimental data. The effect of 
annealing is analyzed on the material constants. 
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1 Introduction 



This study deals with the elastoplastic behavior of injection-molded isotactic polypropy- 
lene (iPP) in isothermal uniaxial tests with small strains. Plastic deformations and yield- 
ing of iPP have been the focus of attention in the past decade, see, e.g., |], |3], £|, ||, || 0, 
H], H , to mention a few. This may be explained by numerous applications of polypropylene 
in industry (ranged from oriented films for packaging to nonwoven fabrics and reinforcing 
fibres) . 

Isotactic polypropylene is a semi-crystalline polymer containing three different crys- 
tallographic forms ||1Q|| : monoclinic a crystallites, hexagonal (3 structures, orthorhombic 
7 polymorphs, and "smectic" mesophase (arrays of chains with a better order in the 
longitudinal than in transverse chain direction). At rapid cooling of the melt (which is 
typical of the injection-molding process), a crystallites and smectic mesophase are mainly 
developed, whereas (3 and 7 polymorphs are observed as minority components JTTJ . 



The characteristic size of a spherulites in injection-molded specimens is estimated as 
100 to 200 /im |TTJ. These spherulites consist of crystalline lamellae with thickness of 10 
to 20 nm [[| [Tj|. A unique feature of a spherulites in iPP is the lamellar cross-hatching: 
development of transverse lamellae oriented in the direction perpendicular to the direction 
of radial lamellae in spherulites [fjl| [12j . 



The amorphous phase is located (i) between spherulites, (ii) in "liquid pockets" JL3 
between lamellar stacks inside spherulites, and (iii) between lamellae in lamellar stacks. 
It consists of (i) mobile chains between spherulites, in liquid pockets and between radial 
lamellae inside lamellar stacks, and (ii) chains with restricted mobility (the so-called "rigid 
amorphous fraction" Jl3|) in regions bounded by radial and tangential lamellae. 

Stretching of iPP specimens induces inter-lamellar separation, rotation and twist of 
lamellae, fine and coarse slip of lamellar blocks and their fragmentation |5|, chain 
slip through the crystals, sliding and separation of tie chains [§], |7|], and activation of 
the rigid amorphous fraction induced by disintegration of transverse lamellae. At large 
strains, these transformations result in cavitation, formation of fibrills and stress-induced 
crystallization fll4 



Annealing and isothermal crystallization of iPP have attracted an essential attention 
in the past five years, see, e.g., jR], [HJ, [Uj, |16|, [T7|, [HJ, Dramatic changes are ob- 
served in DSC (differential scanning calorimetry) traces of isotactic polypropylene driven 
by annealing at elevated temperatures (in the range from 110 to 170 °C). It is found that 
annealing in the low-temperature interval (between 110 and 150 °C) results in (i) a mono- 
tonical increase in the melting peak |10| and (ii) formation of a broad low-temperature 
shoulder (second endotherm) on a melting curve; the intensity of this shoulder grows with 
annealing temperature |18| . Annealing in the high-temperature interval (between 150 and 
170 °C) causes transformation of the second endotherm into the main peak [|I8|| , which 



may be attributed to a second-order phase transition in the crystalline phase f!9 l. The 
critical temperature corresponding to this transition lies in the region between 157 W\ 
and 159 °C ||19|| . The morphological analysis reveals that transformations of the melt- 
ing curves caused by annealing are accompanied by changes in the crystalline structure 
(a pronounced reduction in the level of cross-hatching with an increase in the annealing 
temperature). 
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Although the effect of annealing on the micro-structure of crystallites has been studied 
in detail, a little is known about the influence of thermal treatment on the mechanical be- 
havior of iPP. Our studies on the viscoelastic response of isotactic polypropylene annealed 
in the low-temperature interval (between 110 and 130 °C) revealed that the relaxation 
process is strongly affected by annealing p0| . The aim of the present work is to evaluate 
the effect of thermal treatment on the elastoplastic response of iPP. 

An important shortcoming of conventional elastoplasticity theories for solid polymers, 
see, e.g., [^T], |22J, ^3], [21, ^3J, is that they are grounded on the concept of yield surface. 



According to that approach, below the yield strain e y (which is associated with the point 
of maximum on a stress-strain curve), a specimen demonstrates the nonlinear elastic 
(or viscoelastic) response. This implies that the loading and unloading paths of stress- 
strain diagrams should coincide (provided that the strain rate is so large that the stress 
relaxation phenomenon may be disregarded). The experimental data presented in Section 
2 show that this assumption is inapplicable to isotactic polypropylene, whose stress-strain 
curves for loading and unloading noticeably differ from each other even when the maximal 
strain in a cyclic test is far below the yield strain. 

The objective of this paper is to develop constitutive equations for the elastoplastic 
behavior of a semicrystalline polymer that do not employ the concept of yield surface and 
to apply these relations in fit observations in tensile cyclic tests with maximal strains in 
the sub-yield interval of deformations. 

To make a constitutive model tractable from the mathematical standpoint, we treat a 
semicrystalline polymer as an ensemble of meso- regions (MRs). A meso-domain is thought 
of as an equivalent network of macromolecules bridged by junctions. The mechanical be- 
havior of a MR is associated with that of the amorphous phase, whereas the links between 
meso-domains that transmit the macro-strain to individual MRs reflect the responses of 
crystallites. 

A meso-region is treated as a linear elastic medium, whereas deviations of the stress- 
strain curves in tensile tests from straight lines (that describe the response of a linear 
elastic solid) are attributed to sliding of junctions between chains with respect to their 
reference positions in a stress-free material. This is an essential simplification of the 
mechanism of micro- deformations, which is tantamount to the assumption that only the 
amorphous phase is altered under active loading, whereas the crystalline phase remains 
unchanged. Elastic deformation of spherulites, inter-lamellar separation and fine slip of 
lamellar blocks are not introduced into the model explicitly, but are taken into account 
implicitly in terms of "average" parameters that characterize sliding of junctions in meso- 
domains. This "generalized" sliding process is described by a plastic strain, e p i, whose 
rate of growth is assumed to be proportional to the rate of straining. 

At unloading, junctions between chains in MRs move back to their initial positions 
with a decrease in the macro-strain e. Because deformation of a semicrystalline polymer 
results in micro-fracture of crystallites, the plastic strain e p i (that decreases together 
with the macro-strain e) is not sufficient to describe morphological transformations under 
unloading. To account for coarse slip and disintegration of lamellar blocks at unloading, 
another plastic strain, e P 2, is introduced, which increases in time with a rate proportional 
to the rate of work of external forces. 

A similarity may be noted between splitting the plastic strain, e p , into two components, 
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e p i and e P 2, where the latter quantity is altered at unloading only, and a phenomenolog- 
ical model for the Mullins effect in particle- reinforced elastomers p6 |. According to the 



Ogden-Roxburgh concept, the difference between the stress-strain curves at active de- 
formation and unloading is attributed to some damage parameter that changes along the 
unloading path of a deformation history only. An advantage of our model compared to 
the Ogden-Roxburgh approach is that the plastic strain, e p2 , obeys a conventional flow 
rule in elastoplasticity, whereas the damage parameter introduced in is governed by 
a kinetic equation whose physical meaning is not clear for semicrystalline polymers. 

The exposition is organized as follows. Observations in uniaxial loading-unloading 
tensile tests are reported in Section 2. Kinetic equations for the plastic strains, e p i and 
e P 2, are introduced in Section 3. Constitutive equations for a semicrystalline polymer 
at isothermal uniaxial loading are derived in Section 4. Adjustable parameters in the 
stress-strain relations are determined in Section 5 by matching the experimental data. A 
brief discussion of our findings is presented in Section 6. Some concluding remarks are 
formulated in Section 7. 



2 Experimental procedure 

Isotactic polypropylene (Novolen 1100L) was supplied by BASF (Targor). ASTM dumb- 
bell specimens were injection molded with length 148 mm, width 10 mm and thickness 
3.8 mm. Uniaxial tensile tests were performed at room temperature on a testing machine 
Instron-5568 equipped with electro-mechanical sensors for the control of longitudinal 
strains in the active zone of samples. The tensile force was measured by a standard load 
cell. The longitudinal stress, a, was determined as the ratio of the axial force to the 
cross-sectional area of stress-free specimens. 

Four series of experiments were performed. In the first series, the specimens were used 
as produced. In the other series, the samples were annealed prior to testing for 24 h at 
the temperatures 120, 140 and 160 °C and slowly cooled by air. To minimize the effect of 
physical aging, mechanical tests were carried out at least one day after thermal treatment. 

In any test, a specimen was stretched with a cross-head speed of 10 mm/min (which 
corresponds to the Hencky strain rate e# = 2.09- 10~ 4 s _1 ) up to the maximal strain, e max , 
and unloaded with the same cross-head speed to the zero stress. The chosen cross-head 
speed ensures nearly isothermal experimental conditions p7 [, on the one hand, and it is 



sufficiently large to disregard the viscoelastic effects, on the other (the maximal duration 
of a loading test does not exceed 1 min). 

Any series of experiments consisted of 6 tests with the maximal strains e max = 0.02, 
0.04, 0.06, 0.08, 0.10 and 0.12. Each test was performed on a new specimen. The interval 
of strains under consideration includes the entire sub-yield region (the yield strain for a 
non-annealed iPP e y = 0.13, according to the supplier). We confine ourselves to relatively 
small deformations, because the purpose of the study is to approximate experimental data 
by using constitutive equations with small strains. 

For the sake of brevity, we present only the stress-strain curves for specimens annealed 
at 140 °C. The engineering stresses, a, is plotted versus the engineering strain, e, in 
Figures 1 to 6. The stress-strain diagrams for non-annealed specimens, as well as for 
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samples annealed at 120 and 160 °C have a similar shape. The following features of the 
experimental curves are worth to be mentioned: 



1. at all maximal strains, e max , the stress-strain curves corresponding to the loading 
and unloading paths substantially differ from each other, 

2. the unloading curves are strongly nonlinear, 

3. the residual strain (measured at the instant when the stress vanishes) noticeably 
increases with the maximal strain, e max . 

3 A micro- mechanical model 

A semicrystalline polymer is treated as an equivalent network of chains. The network is 
assumed to be strongly heterogeneous. This heterogeneity is attributed to an inhomogene- 
ity of interactions between chains in the amorphous phase and crystalline lamellae with 
various lengths and thicknesses. The network is modelled as an ensemble of meso-regions 
(MRs) with arbitrary shapes. The characteristic length of a MR substantially exceeds the 
radius of gyration for a macromolecule, and it is noticeably less than a size of a sample. 
Deformation of a specimen induces two processes in the network: 

1. Sliding of junctions (physical cross-links and entanglements that bridge chains in 
the network) with respect to their positions in the stress-free state. 

2. Sliding of meso-domains in the ensemble with respect to each other. 

Sliding of junctions describes non-affine deformation of the network. This process is 
determined by a plastic strain e p \. Sliding of MRs with respect to each other is described 
by a plastic strain e P 2. In a semicrystalline polymer, the strain e p i reflects sliding of 
junctions in the amorphous phase, slippage of tie chains and fine slip of lamellar blocks. 
The strain e p2 characterizes coarse slip of lamellar blocks and their disintegration. 

The total plastic strain, e p , equals the sum of the plastic strains driven by sliding of 
nodes and mutual displacements of meso-domains, 



We suppose that meso-domains are connected by links that transmit the macro-strain, 
e, to individual MRs. This implies the conventional hypothesis that the macro-strain, e, 
equals the sum of the elastic strain in meso-regions, e e , and the plastic strain, e p , 



To simplify the analysis, we assume that the elastic strains in various MRs coincide. 

Deformation of a specimen results in evolution of the plastic strain, e p i, both at the 
stages of active loading and unloading. We suppose that the rate of changes in the plastic 
strain, e p i, is proportional to the rate of changes in the macro-strain, e, 



e p — e pl + e p2- 



(1) 



e = e c + e p . 



(2) 



de 



pi 



(t) = tp(e e (t))-(t), 



(3) 



d,t 
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where the coefficient of proportionality, if, depends on the elastic strain e e . It is as- 
sumed that the function ip{t c ) vanishes at the zero elastic strain, <f(0) = 0, monotonically 
increases with e e , and tends to some constant a e [0, 1] for relatively large elastic strains, 

lim if(e e ) = a, 

where a is the rate of sliding of junctions in MRs corresponding to a developed plastic 
flow. The inequality a > means that junctions slide in the direction that is determined 
by the macro-strain e. The condition a < 1 ensures that the rate of sliding does not 
exceed the rate of straining. 

To approximate experimental data, we employ the function 



<f(e e ) = a 1 - exp(- . (4) 



This function is determined by two adjustable parameters, a and e, where the quantity e 
characterizes how "large" an elastic strain, e e , is. 

With reference to |], [|], we suppose that under active loading in the sub-yield region, 
lamellar fragmentation does not occur and the plastic strain e P 2 vanishes, 

= (|W>°> ^)> )- ( 5 ) 

According to Eq. (5), the entire dissipation of energy at active deformation is attributed 
to sliding of junctions with respect to their reference positions, whereas displacements of 
MRs with respect to each other are disregarded. 

It is assumed that under unloading, meso-domains slide with respect to each other as 
they are driven by a positive macro-stress o. These mutual displacements of MRs are 
characterized by a plastic strain, e P 2, that grows with time t. An increase in e P 2 during 
unloading reflects coarse slip and fragmentation of deformed lamellae. The evolution of 
the plastic strain, e P 2, is described by the flow rule 

Equation (6) means that the rate of changes in e P 2 is proportional to the rate of work of 
external loads 

At) = -*(t)j t (t), 

where the sign "— " accounts for the opposite directions of the longitudinal stress and the 
strain increment. With reference to | 28| , we suppose that the coefficient of proportionality, 



K > 0, is a function of the maximal plastic strain, e pi , reached: K = K(e° 1 ), where 

epiO) = maxe p i(r). 
To match observations, we use the linear function 

K = K + K 1 e° pl , (7) 
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where K m (m = 0, 1) are adjustable parameters. Equation (7) means that the larger 
is the plastic strain, e p i, driven by fine slip of lamellar blocks under active loading, the 
higher is the rate of the plastic strain, e P 2, that reflects coarse slip and fragmentation of 
lamellae at unloading. 

Separation of tie chains from lamellae and disintegration of lamellar blocks result 
in a decrease in the number of MRs to which the macro-strain, e, is transmitted by 
surrounding meso-domains. This decrease is attributed to two different processes: (i) 
mechanically-induced separation of individual MRs from the ensemble, and (ii) screening 
of meso-domains by stacks of disintegrated lamellae. 

To describe evolution of an ensemble of meso-domains, we introduce the average num- 
ber of MRs, No, per unit mass of a virgin specimen (where all MRs are connected with 
one another) and the average number of MRs, N(t), in the deformed specimen at time 
t > (where some meso- regions are separated from the ensemble). 

At the stage of active loading in the sub-yield region, separation of MRs from the 
ensemble does not take place, which implies that N(t) remains constant, 

^(t)=0 (^(f)>0, a(t)>6), N(0)=N . (8) 

Changes in the function N(t) at unloading are governed by a first order kinetic equation, 
according to which the relative number of MRs separated from the ensemble per unit time 
is proportional to the increment of the plastic strain, e P 2, 

-m>= K > (§««»■ <*>*»)■ < 9 > 

By analogy with Eq. (6), the coefficient k > is treated as a function of the maximal 
plastic strain e° x : k = «(ep\). The pre-factor k in Eq. (9) is assumed to be rather large for 
small plastic strains, e p i, i.e., when straining of a specimen produces no fine slip of lamellar 
blocks, and it monotonically decreases with e p i. This decrease is attributed to alignment 
of lamellar blocks driven by their fine slip, which, in turn, diminishes the probability of 
formation of disordered lamellar stacks at unloading. These stacks do not transmit the 
macro-strain, e, to the amorphous domains surrounded by them, which results in isolation 
of these regions from the ensemble (in a way silimar to the formation of regions of occluded 
rubber in particle-reinforced elastomers [f29f ). 

To fit experimental data, we adopt the function 

e ° 

K = K + K 1 expy — —J, (10) 

where K m (m = 0, 1) and e* are adjustable parameters. Equation (10) implies that k 
monotonically decreases from the initial value k,q + «i and tends to the limiting value kq 
at relatively large plastic strains, e p i. The quantity e* indicates how "large" the maximal 
plastic strain, e pl , is. 
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4 Constitutive equations 



Under isothermal uniaxial deformation, a meso-domain is treated as a linear elastic solid 
with the mechanical energy 



where the constant fx > is the average rigidity of a MR. Neglecting the energy of 
interaction between meso-regions, we calculate the strain energy density per unit mass of 
a polymer as the sum of the mechanical energies of MRs, 



W{t) = ^N{t)el 



It follows from Eqs. (1) to (3) and (11) that the derivative of the function W with respect 
to time is given by 



dW , , fxdN 
it) 



dt y ' 2 dt 
The Clausius-Duhem inequality reads 



(t)el(t) +fxN(t)e c (t)[(l - - 



dt 



;i2) 



, s dW , . 1 , ,de , , 

Q{t) = — (i) + r (t) s (t) £ °- 

where p is mass density, and Q is internal dissipation per unit mass. Substition of Eq. 
(12) into this formula results in 



Q(t) 



1 T a(t) - P fxN(t)e e (t)(l - <p(e e (t)))} ^(t) + Y(t) > 0, 



where 



rde 



Y(t)=fxN(t)e e (t) =£(t) 



dt 

1 dN 



(t)e e (t) 



dt w 2N(t) dt 

Equating the expression in square brackets in Eq. (13) to zero, we find that 



a(t) = E(t)e e (t) l-y?(e c (t)) 



with 



E(t) = pfiN(t). 

Equations (5), (8) and (14) imply that the function Y vanishes at active loading, 

Y(t)=0 (j t (t)>0, a(t)>0). 
It follows from Eqs. (6), (9), and (14) to (16) that at unloading, 



(13) 
(14) 

(15) 
(16) 

(17) 



a 2 (t) de 



de 



<p(te(t)) dt 



(18) 
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Combining Eqs. (13), (15) and (16), we find that the internal dissipation per unit mass, 
Q(t), coincides with Y(t). According to Eqs. (17) and (18), the function Y(t) is non- 
negative. This means that the Clausius-Duhem inequality is satisfied for an arbitrary 
deformation program, provided that the stress, a, is given by Eq. (15). 

It follows from Eqs. (8), (9) and (16) that the elastic modulus, E, obeys the differential 
equations 

^§(t) = 0, E(0) = E (§(<)><>, »(*)>0 

^>-4f« (><°- '<«>*<>)■ < 19 > 

where Eq = fipNo is the elastic modulus for a virgin specimen. 

The assumption that elastic moduli are strongly affected by mechanical factors is 
widely used in elastoplasticity theories for geomaterials [0, |3l| . For granular media, 
this effect is attributed to the growth of voids between particles under active loading 
and to the contractive response at unloading (which means that elastic moduli become 
functions of the volumetric strain In the present model, the plastic strain, e P 2, 

plays a role similar to that the mean elastic strain plays in elastoplasticity theories for 
granular materials: it characterizes the level of disorder produced by mechanically-induced 
lamellar fragmentation. The difference between our approach and previous studies is that 
for an arbitrary time-dependent loading program, the elastic modulus, E, depends not 
on the current strain, but on the entire history of deformation [because Eqs. (19) cannot 
be integrated explicitly when k is a function of the maximal plastic strain, e° 1; reached 
under active loading]. 

Uniaxial deformation of a semi- crystalline polymer is described by Eqs. (1), (2), (4) to 
(6), (15) and (19). Any stress-strain curve for cyclic loading is determined by 5 adjustable 
parameters: 

1. the initial elastic modulus Eq, 

2. the rate of a developed plastic flow a, 

3. the strain e that characterizes transition to a steady plastic flow, 

4. the rate of sliding of MRs with respect to each other K = K^e^), 

5. the rate of separation of meso-domains from an ensemble k = ^(e^). 

This number is quite comparable with the number of material constants in other consti- 



tutive models in elastoplasticity, see, e.g., [22, 23, 25]. It should be noted, however, that 



most of these concepts fail to adequately describe the mechanical behavior of polymers in 
cyclic tests |32] . 



An important advantage of the stress-strain relations (1), (2), (4) to (6), (15) and (19) 
is that 3 constants, Eq, a and e, are found by fitting experimental data for the loading 
path of a stress-strain curve, whereas the other two parameters, K and k, are determined 
by matching observations for the unloading path. 
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The dependences of K and k on the maximal plastic strain, e° pl , are described by Eqs. 
(7) and (10). To find the coefficients K m , K m and e* is these equations, the quantities 
K and k (determined by matching observations in a series of cyclic tests with various 
maximal strains e max ) are approximated by functions (7) and (10). 



5 Fitting of observations 

We begin with the approximation of the stress-strain diagrams for active loading. It 
follows from Eqs. (1), (2), (4), (5), (15) and (19) that the stress, a, is given by 

<r(e) = E (e - e p i){l - a[l - exp(-^=^)] }, (20) 
where the plastic strain, e p i, satisfies the nonlinear differential equation 

^(6) = a[l-exp(-^i)], e pl (0) = 0. (21) 



The loading path of any stress-strain curve is determined by 3 material constants: E , 
a and e. To find these quantities, we fix some intervals [0, a max ] and [0, £ max ], where the 
"best-fit" parameters a and e are assumed to be located, and divide these intervals into J 
subintervals by the points a* = iAa and Ej = jAe (i, j = 1, . . . , J) with Aa = a ma , x /J and 
Ae = £ max /J. For any pair, {ai,Ej}, Eq. (21) is integrated numerically by the Runge- 
Kutta method with the step Ae = 1.0 • 10~ 5 . Given a pair, {a^e,-}, the elastic modulus, 
E = E (i, j), is found by the least-squares method from the condition of minimum of the 
function 

F(hj) = ^[°"cxp(e m ) - 0"num(e m ) 



2 



where the sum is calculated over all experimental points, e m , on a loading path, a exp 
is the stress measured in a tensile test, and <7 num is given by Eq. (20). The "best-fit" 
parameters a and e are determined from the condition of minimum of the function F on 
the set {ai,Ej — 1, . . . , J)}. 

The material constants Eq, a and e that minimize the discrepancies between the 
experimental data and the results of numerical analysis are found for any stress-strain 
curve independently. Afterwards, we calculate the average values of these quantities 
and their standard deviations for 4 series of tests. These values are presented in Table 
1 for non-annealed specimens and for samples annealed at 120, 140 and 160 °C. The 
table shows that the elastic modulus, Eq, and the rate of developed plastic flow, a, are 
determined with a high level of accuracy (the relative deviations are less than 8 and 13 
%, respectively), whereas the accuracy in determining the strain, e, that characterizes 
transition to a steady-state plastic flow, is rather low (about 29 %). 

To find the quantities K and k, we approximate the unloading paths of the stress- 
strain curves. It follows from Eqs. (1), (2), (4), (6), (15) and (18) that the stress, a, is 
given by 

a(e) = E(e p2 )(e - e pl - e p2 ){l - a[l - exp(~ e ~ Cpl ~ Cp2 )] }, (22) 
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where the elastic modulus, E, reads 



E(e p2 ) = E exp(-Ke p2 )- 
The plastic strains, e p i and e P 2, obey the differential equations 



de 



pi 



de 



cxp 



(- 



-pi 



-p2 



P 2 



e) = -tf(7(c) 



(23) 



(24) 



with the "initial" conditions 



Cpl(Cmax) — Cpl? Cp2(Cmax) — 0. 

To approximate experimental data in a cyclic test with a maximal tensile strain e max , 
we apply an algorithm similar to that used to match the stress-strain curves at active 
loading. We fix some intervals [0, K max \ and [0, /t max ], where the "best-fit" parameters K 
and k are assumed to be located, and divide these intervals into J subintervals by the 
points Ki = iAK and Kj = jAk = 1, . . . , J) with AK = K max / J and Ak = K max / J. 
For any pair, {Ki, Kj}, Eqs. (20) to (22) are integrated numerically by the Runge-Kutta 
method with the step Ae = 1.0- 10~ 5 and with the value of Eq found by fitting the loading 
paths of the stress-strain curves. The "best-fit" parameters K and k are determined from 
the condition of minimum of the function F on the set {Ki, Kj (i, j = 1, . . . , J)}. 

The parameter K is plotted versus e° x in Figure 7. The maximal plastic strain, e° l5 is 
determined by numerical integration of Eq. (21) from e = to e = e max . The experimental 
data are fitted by Eq. (7), where the coefficients K m (m — 0,1) are found by the least- 
squares technique. 

The quantity k is depicted as a function of e° x in Figure 8 together with its approxi- 
mation by Eq. (10). The strain e* is determined by the steepest-descent method, whereas 
the coefficients n m (m = 0, 1) are calculated by the least-squares algorithm. Unlike Figure 
7, where the only curve provides an acceptable approximation of all experimental data, 
adjustable parameters in Eq. (10) are found independently for specimens not-subjected 
to thermal treatment and annealed at temperatures in the low-temperature region, on the 
one hand, and for samples annealed at 160 °C, on the other. 



6 Discussion 

Figure 1 to 6 demonstrate fair agreement between the observations in loading-unloading 
tests with various maximal strains, e max , and the results of numerical simulation. The 
same quality of fitting has been reached also for the stress-strain curves for specimens 
not subjected to thermal treatment and those annealed at 120 and 140 °C. These figures 
confirm that the model correctly describes the mechanical response of semicrystalline 
polymers in cyclic tests. 

Table 1 shows that annealing of iPP in the low-temperature region does not affect the 
elastic modulus, E , and the rate of developed plastic flow, a. Annealing at 160 °C causes 
a noticeable increase (from 20 to 30 %) in the elastic modulus, E , and a growth of the 
rate of steady plastic flow, a, that reaches its ultimate value. Although the increase in 
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a is not rather high (about 6 to 14 %), the fact that the rate of developed plastic flow 
reaches unity means that at relatively large elongations, the elastic strain, e e , does not 
grow, and an increase in the macro-strain, e, is totally compensated by the same increase 
in the plastic strain, e p \. On the contrary, the rate of developed plastic flow in specimens 
not subjected to thermal treatment, as well as in those annealed in the low-temperature 
region is noticeably less than unity, which implies a strong increase in the elastic strain 
under stretching. This difference becomes important for the analysis of stress-strain 
curves at large elongations (appropriate data are not presented): necking of specimens 
not subjected to annealing occurs at the Hencky strains about 0.2, whereas necking of 
specimens annealed in the high-temperature region is not observed at the Hencky strains 
up to 0.6 . 

The growth of the elastic modulus, E Q , of specimens annealed at 160 °C may be at- 
tributed to an increase in perfectness of crystalline lamellae (associated with transition 
from a statistically disordered a.\ phase into an ordered a 2 phase characterized by regu- 
larity in the up and down positions of methyl groups along the chains [|T^J). An increase 
in the rate of developed plastic flow, a, may be ascribed to a decrease in the level of 
cross-hatching in a spherulites with the growth of annealing temperature |T^] and the 
total disappearance of transverse lamellae above the critical temperature T c = 159 °C 

G3- 

Table 1 shows that the parameter e is independent of the annealing temperature 
(within the range of experimental uncertainities). This conclusion appears to be quite 
natural, because e characterizes sliding of junctions between chains in MRs. This quantity 
describes the response of the amorphous phase that is not affected by thermal treatment. 
The value of e (in the range between 3 and 4 %) found by fitting the stress-strain curves 
appears to be rather close to a critical strain, e c , determined by matching relaxation curves 
as the strain at which the rate of relaxation becomes independent of mechanical 



factors. 

Figure 7 reveals that the parameter K is not affected by annealing temperature. At 
first glance, this result seems rather surprising, because K is responsible for mutual dis- 
placements of MRs driven by coarse slip and fragmentation of lamellar blocks. Two ex- 
planations may be provided for our findings. The first is based on the Nitta-Takayanagi 
Q and Meyer-Pruitt [|34] concepts that attribute plastic deformation of semicrystalline 
polymers to slippage of tie chains along lamellae and pulling out of chains from disinte- 
grated lamellar blocks (both processes are associated with transformations in the amor- 
phous phase whose state is not influenced by thermal treatment). The other explanation 
is based on the assumption that an increase in perfectness of crystallites at annealing 
(which makes lamellae stronger) is compensated by disappearance of tangential lamellae 
(which results in weakening of spherulites). Because these morphological transformations 
in iPP lead to opposite changes in the mechanical response, their combined effect may be 
negligible. 

According to Figure 8, the parameter k is not affected by annealing in the low- 
temperature region, but is strongly influenced by high-temperature annealing. The limit- 
ing value of this quantity, kq, that corresponds to a developed plastic flow at unloading, in- 
creases approximately by twice after annealing at 160 °C. Because k characterizes changes 
in the elastic modulus driven by lamellar fragmentation, see Eq. (19), we associate this 
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growth with weakening of spherulites induced by disappearance of cross-hatching. 

The strain, e*, that determines transition to a steady plastic flow at unloading (for all 
specimens, but those annealed at 160 °C) is very close to the strain, e, that characterizes 
transition to a developed plastic flow at active loading. This implies an opportunity to 
reduce the number of adjustable parameters in the constitutive equations by replacing 
e* by e. However, an additional analysis of the physical basis for this simplification is 
required. 

7 Concluding remarks 

Four series of tensile loading-unloading tests have been performed on injection-molded 
isotactic polypropylene at room temperature. In the first series, specimens were used as 
produced, whereas in the other series, samples were annealed for 24 h at 120, 140 and 160 
°C prior to testing. 

Constitutive equations have been derived for the elastoplastic response of semicrys- 
talline polymers. A polymer is treated as an equivalent network of chains bridged by 
junctions. The network is assumed to be strongly heterogeneous, and it is modelled as an 
ensemble of meso-regions linked with each other. Under active loading in the sub-yield 
region of deformations, junctions between chains in MRs slide with respect to their refer- 
ence positions (which describes sliding of tie chains and fine slip of lamellar blocks). At 
unloading, this non-afrme deformation of the network is accompanied by displacements of 
meso-domains with respect to each other (which describes coarse slip and fragmentation 
of lamellar blocks). Destructure of lamellae results in separation of some MRs from the 
ensemble (driven by breakage of links between isolated meso-domains and the network and 
screening of macro-strain in these meso-domains by surrounding stacks of disintegrated 
lamellae. 

Kinetic equations are proposed for the rates of two plastic strains that describe these 
processes and for the rate of separation of MRs from an ensemble. Stress-strain relations 
for isothermal uniaxial deformation are developed by using the laws of thermodynamics. 
The constitutive equations are determined by 5 adjustable parameters that are found by 
fitting the experimental data. Fair agreement is demonstrated between the experimental 
stress-strain curves and the results of numerical simulation. 

The following conclusions are drawn: 

1. Annealing in the low-temperature region does not affect the material constants that 
reflect the elastoplastic response of iPP. 

2. Annealing in the high-temperature region results in an increase in the elastic mod- 
ulus E (which is attributed to the growth of the perfectness of crystals induced by 

— > oi2 transition). 

3. Annealing in the high-temperature region causes the growth of the rate of devel- 
oped plastic flow a (that reaches its ultimate value) and an increase in the rate of 
separation of meso-regions from an ensemble k (these changes are associated with 
disappearance of transversal lamellae). 
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4. The rate of plastic strain, K, linearly grows with the maximal plastic strain, e° 1; 
which means that sliding of junctions in MRs under active loading activates coarse 
slip and fragmentation of lamellae at unloading. The activation process appears to 
be independent of the perfectness of crystallites. 

5. The strains e and e* that characterize transitions to steady plastic flows at active 
loading and unloading, respectively, are rather close to each other and are weakly 
affected by thermal treatment. This result confirms the hypothesis [|], [34[ that 
elastoplastic deformation of a semicrystalline polymer in a sub-yield region is mainly 
associated with transformations in the amorphous phase. 
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List of figures 



Figure 1: The stress a MPa versus strain e in a tensile loading-unloading test with 
the maximum strain e max = 0.02. Circles: experimental data on a specimen annealed at 
T = 140 °C. Solid line: results of numerical simulation 

Figure 2: The stress a MPa versus strain e in a tensile loading-unloading test with 
the maximum strain e max = 0.04. Circles: experimental data on a specimen annealed at 
T = 140 °C. Solid line: results of numerical simulation 

Figure 3: The stress a MPa versus strain e in a tensile loading-unloading test with 
the maximum strain e max = 0.06. Circles: experimental data on a specimen annealed at 
T = 140 °C. Solid line: results of numerical simulation 

Figure 4: The stress er MPa versus strain e in a tensile loading-unloading test with 
the maximum strain e max = 0.08. Circles: experimental data on a specimen annealed at 
T = 140 °C. Solid line: results of numerical simulation 

Figure 5: The stress er MPa versus strain e in a tensile loading-unloading test with 
the maximum strain e max = 0.10. Circles: experimental data on a specimen annealed at 
T = 140 °C. Solid line: results of numerical simulation 

Figure 6: The stress er MPa versus strain e in a tensile loading-unloading test with 
the maximum strain e max = 0.12. Circles: experimental data on a specimen annealed at 
T = 140 °C. Solid line: results of numerical simulation 

Figure 7: The parameter K MPa -1 versus the maximal plastic strain e° pl . Symbols: 
treatment of observations. Unfilled circles: non-annealed specimens; filled circles: speci- 
mens annealed at T = 120 °C; asterisks: specimens annealed at T = 140 °C; diamonds: 
specimens annealed at T = 160 °C. Solid line: approximation of the experimental data 
by Eq. (7) with K = 0.0066 and K x = 0.2659 

Figure 8: The dimensionless parameter k versus the maximal plastic strain e^. Symbols: 
treatment of observations. Unfilled circles: non-annealed specimens; filled circles: speci- 
mens annealed at T = 120 °C; asterisks: specimens annealed at T = 140 °C; diamonds: 
specimens annealed at T = 160 °C. Solid lines: approximation of the experimental data 
by Eq. (10). Curve 1: k = 57.96, «i = 215.83, e* = 1.23 • 10~ 2 ; curve 2: k = 106.68, 
«i = 1184.71, e* = 3.83- 10~ 3 
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Table 1: Adjustable parameters Eq GPa, a and e for specimens annealed at various 
temperatures T °C (the values in parentheses indicate standard deviations) 



T 


E 


a 


e-10 2 


non-annealed 


1.81 (0.11) 


0.93 (0.04) 


3.84 (0.53) 


120 


1.75 (0.08) 


0.87 (0.09) 


3.52 (0.50) 


140 


1.71 (0.14) 


0.91 (0.12) 


4.18 (1.23) 


160 


2.18 (0.12) 


0.99 (0.01) 


3.84 (0.33) 
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